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ABSTRACT 


A ten-level primitive equation ocean circulation model 
is used to investigate the formation and propagation of baro- 
clinic wavelike responses to annual variations in wind stress 
monreing. At 12.2°N a thermal response to Ekman pumping is 
@ee-aved tO be in phase across the entire longitudinal extent 
meee model Pacific Ocean. A Zonal fluctuation in interior 
ocean current systems is also suggested. Eastern boundary 
disturbances not related to local forcing are observed to 
propagate poleward. At 6.1°N, main thermocline displacements 
are observed to propagate westward as the result of vertical 
advection induced by Ekman pumping. Surface layer tempera- 
ture disturbances also exist which are produced by meridional 
advection of mean temperature (colder to the south) by sur- 
face Ekman currents. Eastern boundary disturbances also 
exist at this latitude and extend to a deeper level. All 
disturbances obServed in this investigation are considered 
to be permanent features of the model circulation. Grid size 
limitations preclude detailed investigation of the observed 


boundary disturbances with this model. 





TABLE OF CONTENTS 


ie LITRODUCTION = 25 10 
MED ATA SOURCES AND PROCESSING - - - - - - - - - - 1 
MPM—ZOSSERVATIONS AND RESULTS - - - -------- - 16 
GOINGS WON 5 Sl 
PCP URENCES=6= 9-2 -9= == = = = -=-=-+--- - 65 


INITIAL DISTRIBUTION LIST - ------+-7+77-7-7-- oe 





Loy 


ie. 


feo 2b ear fouURES 


Meridional GausSian smoothing function used 
over the surface wind stress field as de- 
scribed in text - - - - - - ------7- 7-7-7 


Annual mean surface wind stress vectors over 

the model domain. Stress vectors are plotted 

at aiternate grid points with the length of 

mae arrow directiy proportional to the magni- 
tude of the stresS - ---+--+-+-+--+-+-+- --- 


Sur of the surface wind stress for a typical 
serine month (April) as used in the model. 
Contour interval is 5 x 1073 N/m’. Hatched 

areas represent positive values - - - ------ 


same as Figure 3 except for a typical fall 
Prec Oeremen je = 9= = = — — = = ee a es 


Time-longitude map of the departure from the 
mimual mean of the zonal (u') fiow across 

me 2°N at 162.5 m (level 5). Contour inter- 
Waeets O.4¢ cm/s. Hatched areas represent 
positive values - - ---+-+-+--+-+-+-+-+-+-+-+--+- 


= 


pame as Figure 5 except for meridional (v') 
a es ee 


Same aS Figure 5 except for vertical velocity 
er). Contour interval is 2.0 em/day - - - ~ = - 


Time-latitude map of the departure from the 
annual mean of vertical velocity (w') along 

Maes mocdel Castern boundary (i125°W) at level 5. 
Seutour interval is 30.0 cm/day. Hatched 

areas represent positive values - - ---+---+--- 


Same as Figure 5 except for temperature (T'). 
i Veemwmameet aims Ol.) C.= = = = — - = —- = —~ — = 


Same as Figure 5 except for Ekman pumping 
Zech vyweecOntounr Interval is 4.0 cm/day - - ~ - 


Vertical cross section (x,z) of the departure 
from the annual mean of the zonal flow (u') 
across 12.2°N on August 3 of the last model 

Veae. COntour interval is 0.4 cm/s. Hatched 
areas represent positive values - - ------ - 


34 


20 


oe 


a6 


By, 





12. Same as Figure 11 except for meridional (v') 
flow - - =--=+--+--+-4+e8777-77--+-- ee ean es A 


13. Same as Figure 11 except for vertical velocity 
aye.))6©CONtoOur interval 1s 5.0 em/day - - - - - - 45 


14. Same as Figure 11 except for temperature (T'). 
Paaoemntcrgal 15 O0,18°C —- = =. - - - =. = ~ = = = 46 


15. Time depth map of the departure from the 
maaudal mean of the zonal (u') flow at 
latitude 12.2°N and longitude 162°W. Con- 


meur interval 1s 0.5 cm/s. Hatched areas 
represent positive values - - - ------+--- - 47 


16. Same as Figure 15 except for meridional (v') 


flow ------+-+-7--7-7-77772-~°- --- - - - 48 
17. Same as Figure 15 except for vertical velocity 
(w'). Contour interval is 2.0 cm/day - - - - - - 49 


18. Same as Figure 15 except for temperature (T'). 
Miiaetiumintervalei1s O.1°C = === == —- = — — — — 50 


19. Time-longitude map of the departure from the 
annual mean of the zonal (u') flow across 
oN at level 5. Contour interval is 0.4 


cm/s. Hatched areas represent positive values - dl 


feeeeesame as Figure 19 except for meridional (v') 


flow - ----+-+-+-+-2+e-+-7+7+-77-7-27+024- ~ - - 52 
21. Same as Figure 19 except for vertical velocity 

oy). Contour interval is 2.0 em/day - - - - - - 33 
22. Same as Figure 19 except for temperature (T'). 

Contour interval is 0.1%C - -+---+--+-+-+-+- =~ 54 


23. Same as Figure 19 except for Ekman pumping 
velocity. Contour interval is 4.0 cm/day - - - - 959 


Hee Vertical cross section (x,z) of the departure 
Monet ne annial mean of the zonal (u') flow 
across 6.1°N on August 3 of the last model 


year. Contour interval os 0.4 cm/s. Hatched 

areaS represent poSitive values ~ - - =~ =~ = - - -— 56 
25. Same aS Figure 24 except for meridional (v') 

flow ----+-+----- ------7----- - «57 
26. Same aS Figure 24 except for vertical velocity 


eer connour Interval 1S 3.0 cm/day - = - - - - 58 








27. 


28. 


2. 


10 


Si). 


Die 


Same as Figure 24 except for temperature (T'). 
Peoiieericerva is Ol C= = = == --=-=+=---+-+ - og 


Time depth map of the departure from the 

anoueal mean of the zonal (u') flow at 

latitude 6.1°N and longitude 162°W. Con- 

tour interval is 0.5 cm/s. Hatched areas 

represent positive values - -----+-+-+-+-+-- 60 


Same as Figure 28 except for meridional (v') 
iC ee ee Sa a a ee 61 


foes Of the zonal (Tt,') and meridional (T,,') 
components of wind stress as functions of y 


time at latitude 6.1°N and longitude 162°w - - -62 


Same as Figure 28 except for vertical velocity 
(w'). Contour interval is 2.0 cm/day - ---- - 53 


Same as Figure 28 except for temperature (T'). 
Mentour interval is 0.1°C - - ---------- 64 








ACKNOWLEDGEMENTS 


The author wishes to express his appreciation to 
Dr. Robert L. Haney for his guidance, patient understanding 
and counsel in the preparation of this thesis and to 


Dr. C. P. Chang for reading the thesis and making valuable 


S2mments. 








I. INTRODUCTION 


Recent modeling investigations of the strength and 


@@ermal Structure of Tropical North Pacific Ocean currents 
have provided evidence as to the existence of a baroclinic 
response to annual variations in wind stress. Meyers (1978) 
has observed variations in the depth of the 14°C isotherm 
wiaieh, at 6°N, propagate westward at nearly the speed of 
free, non-dispersive Rossby waves. His model mechanisms at 
6°N are Ekman layer divergence (EkKman pumping) as a response 
to local wind stress, and planetary, geostrophic divergence 
inherent on a rotating sphere. At 10°N Meyers observed 
thermocline depth variations which were nearly in phase 
across the Pacific with the Ekman pumping velocity calcu- 
lated from the local wind stress. White (1977) found an 
agreement in phase only at 150°W, and suggested a baroclinic 
response to annual forcing in the form of a long wave which 
results from the superposition of baroclinic Rossby waves 
emanating from the eastern boundary, and the local forced 
response. This resulting long wave propagated at twice the 
Speed of free Rossby waves. 

The purpose of this study was to investigate further 
any low frequency responses to wind forcing in the North 
Equatorial Pacific and to consider whether any such responses 


are consistent with theory. The numerical model used is that 


of Haney (1974) with subsequent improvements involving 








Paeameterization of convective mixing and surface wind stir- 
Memon ( Haney and Davies, 1976) and nonlinear eddy viscosity 
(Haney and Wright, 1975). This model is based on the hydro- 
static and Boussinesq approximations and is a 10-level prim- 
itive equation model in a closed rectangular basin with a 
uniform depth of 4 kilometers. The model basin extends from 
the equator to 65°N and westward for 90° of longitude from 
125°W. Horizontal resolution is approximately 3° of longi- 
midewand 2° of latitude with 33 uniformly spaced ese ict jeve)abaynes 
in both the east-west and north-south directions. The ten 
memeeal levels are placed at 10, 32.5, 62.5, 102.5, 162.5, 
462.5, 900, 1700, and 3100 meters below the undisturbed sea 
Surface. The model has been developed and improved in order 
to study large-scale thermal anomalies observed in the 
Central North Pacific Ocean by the North Pacific Experiment 
(NORPAX). 

The model was "Spun up" for a number of years of 
Simulated time sufficient to allow circulations and thermal 
patterns to establish themselves and then detailed analysis 
of the transient response in the tropical regions was made. 
The wind stress prescribed as a forcing function in the model 
from the equator to 30°N was that used by Meyers (1978) and 
described in detail by Wyrtki and Meyers (1975). North of 
30°N wind stress was calculated from geostrophic wind data 
Obtained from an NCAR data tape of the Northern Hemisphere 


meen avological atlases of Jenne et al., (1969a, b). The 


major region of interest was south of 20°N. This location 











was where Meyers and White observed wavelike disturbances 
and it was reasonable to expect Similar waves in this study. 
By analyzing any observed waves as to their phase relation- 
Ship with the local wind stress curl, their phase speed, 
and their wavelength, possible mechanisms for their genera- 
tion and propagation may be suggested and compared with 
theory. 

By better understanding how the oceanic circulation 
Changes as a result of variant wind forcing, more compre- 


hensive studies of large-scale air-sea interaction and 


Glimatic changes may be undertaken. 








II. DATA SOURCES AND PROCESSING 


mace trom two different sources, one covering the region 
Memeo of 30°N, and one covering the region south of 30°N , 
were combined to produce a realistic wind stress field which 
was used as a forcing function in this modeling investiga- 
tion. The first data base was that processed and explained 
in detail by Wyrtki and Meyers (1975). This monthly aver- 
M@edewind stress data, sorted into 2° latitude by 10° longi- 
tude fields encompassing 29°S to 29°N and 125°E to 75°W, 
was generously provided by Dr. Steve Piascek of NORDA. Only 
the data which was appropriate for use by the model, i.e. 
Meme data from the equator northward and from 145°E to 125°W, 
were used in this study. Wyrtki and Meyers calculated these 


stress data for each. observation according to the formula 
(ee, 1°} = 


where t iS the wind stress, rectangular component positive 


to the east (A) and north (6), ove FS themeonstant SUurrtace air 


Semoity of 1.2 kg/m? , Cq is the drag coefficient assumed con- 


: d 
; Ua Buc) 7 eine 
| +, 9 2.3 
the surface wind components, and pve = (u, + We Ve a he 


stant in the trade wind zone at 1.5 x 10. 


magnitude of the surface wind. A quadratic interpolation 
scheme was used to transfer the Wyrtki-Meyers monthly wind 


stress component arrays, extending from the equator to 29ON, 


to the model gridpoints. 








The second source of data was an NCAR tape of monthly 
mean climatological surface geostrophic wind component data 
meme Ct al., 1969a, 5). These wind components had been 
sorted into 5° latitude by SS longitude fields encompassing 
the model latitudes and ranging in longitude from 120°R to 
80°W. These data were used to calculate surface wind stress 
using the same formula as above. However, an atmospheric 


meee coefficient of 1.3 x 107° 


Wa Soc CmeBUot igen eet all’, 
1971) and a 10% reduction in magnitude from geostrophic was 
assumed. The resulting stress components for the model 
domain north of 30°N were interpolated to model gridpoints. 

By combining the stress component fields computed and/or 
interpolated from the two data sources, 33 x 33 fields en- 
compassSing the model's geographical domain for individual 
months were obtained. The Wyrtki-Meyers stress data, in- 
cluded from the equator to DEN. 1S based on actual wind 
observations in the tropics, while the NCAR data, from 30°N 
to the northern boundary, is based on geostrophic wind cal- 
culations from monthly mean sea-level pressure data. It is 
felt that the resulting stress field represents the best 
total stress pattern that is presently available from monthly 
mean data over the North Pacific Ocean. 

Although the major feature of the two stress fields 
matched quite well at the common latitude of BON. a smooth-— 
ing technique was employed to remove small scale irregular- 
ities and improve horizontal consistency. A weighted three- 


point meridional average given by 


14 








*« 


was selected with the weight, ee Becausstam fuMmetion Of the 
Peemedde (1 1S the longitudinal gridpoint index and j is the 
iaetcudinal index). is detailed and illustrated in 
Figure 1. At the interface between the two data sources, 
30°N, W, approaches one, and at the latitudinal boundaries 
of the model it approaches zero. Essentially, then, at the 
interface the smoothing process was a 4, 4, 4% meridional 
filter, and near the boundaries was negligible. The merid- 
lonal smoothing was performed twice with boundary values 
being unchanged. Additionally, a #, #, # monthly time 


filter was used once over the entire fields and is given by 


S = f 
ee ke) * 874k * 2,5, n+1)> 


where k indicates the month. 

The ocean circulation model (Haney Cea toe) was 
integrated over a total of nine complete annual cycles start- 
ing from the results of the previous model. The only change 
from the previous model is the introduction of the new wind 
Stress fields as described above. After nine years of inte- 
gration it was clear that the upper layers of the ocean 
model, i.e. those most strongly affected by the annual 
Variation in wind forcing, were essentially repeating them- 
selves one year after another. The results analyzed below 
are those from the eighth and ninth years of the nine year 


integration. Primary attention is devoted to the annual 


variation of model variables, less to the annual mean. 








III. OBSERVATIONS AND RESULTS 


The trade wind field over the Pacific Ocean is one of 
the largest and most consistent wind fields on earth although 
little is actually Known about its time and Spacial fluctua- 
tions (Wyrtki and Meyers, 1976). Wind-driven ocean circula- 
tion studies are largely dependent on wind stress field 
Mmeaeum@actlions Since the curl of the wind stress, and thus 
Ekman surface layer divergence or pumping, is a very crucial 
forcing parameter. Accordingly, maps of the surface wind 
Stress and of its curl were analyzed for Significant fluc- 
tuations which may be related to observable oceanic 
Seeocurbances. 

Figure 2 is a depiction of the annual mean wind stress 
vectors over the model domain. The region of primary inter- 
fein this investigation is that south of 20°N. As Shown , 
the Intertropical Convergence Zone (ITCZ) has an annual 
mean position at about gn. Wyrtki and Meyers (1975) have 
Shown that the curl of the wind stress and thus the Ekman 
pumping velocity is relatively insensitive to changes in wind 
Strength and depend mostly on the locations of strong merid- 
mogat shear, which show latitudinal variation throughout the 
year. In the northern hemisphere winter and spring, the 
northeast trades are at their strongest and the southeast 


trades at their weakest. The ITCZ lies about 3° south of 


itS annual mean poSition and as seen in Figure 3, a map of 








memes tress Curl for a typical spring month, regions of 
strong positive wind stress curl appear between 5°N and 10°N. 
Two zones of particularly strong positive curl are apparent, 
one between 155°E and 175°E and another between 155°W and 
140°W. In the northern hemisphere summer and fall, the 

ITCZ shifts northward to about 13°N and the zone of Strongest 
cyclonic curl is located between 10°N and 15°N as seen in 
Figure 4. 

The strongest meridional shear in the Surface wind stress 
1s observed in a band slightly north of the ITCZ which passes 
Over the annual mean position of the ITCZ twice per year. 
Therefore the amplitude of the first annual variation of the 
Meer is largest at the latitudinal extremes of the ITCZ and 
small at the annual mean position. Meyers (1978) computed 
the amplitude of the first annual harmonic of the EkKman pump- 
ing velocity and the depth of the 14°C isotherm in the North 
Pacific and found that both reached maxima at 6°N and at 10°N. 
Not Surprisingly, these latitudes coincided with the zones of 
strongest cyclonic wind stress curl. AS mentioned previously, 
at 10°N the Ekman pumping velocity and the depth of the 1a) 16 
meonnerm fluctuated in phase across the Pacific while at 6°N 
the temperature fluctuations propagated westward from a 
meeron Of resonant forcing in the Eastern Pacific. In order 
to analyze results which could best be compared with those 
of Meyers (1978), model latitudes of 6.1°N and 12.2°N were 


chosen for study in this investigation. 


lias, 








a. Results at 12.2°N 

Meyers' (1978) result at LO°N waS an oceanic tempera- 
ture response which was nearly in phase (a two month zonal 
lag) with the local Ekman pumping velocity across the entire 
Pacific Ocean. White (1977) theorized that the annual fluc- 
Bien at any interior ocean location from 10°N to 20°N 
depends on the fluctuating wind field over the entire east- 
ward ocean with the phase speed and amplitude strongly related 
to distance from the eastern boundary. In a bounded ocean, 
the presence of the eastern boundary can Significantly alter 
the phase distribution of the oceanic response to Ekman pump- 
ing by acting as a Monesnt nated SORECC Ole bDaroc Linic 
planetary waves which propagate boundary responses westward 
into the interior, drastically altering the phase relation 
Pmeaweted by local forcing only (White, 1977). 

Model level 5, corresponding to a depth of 162.5 
meters, was chosen for detailed study Since it 1S the model 
level which most closely corresponds with the observed depth 
of the 14°C isotherm (Meyers, 1978). This level iS above the 
intermediate layers where the vertical temperature gradient 
1s weak and below the level of the maximum density gradient, 
whose stability shields it from temperature fluctuations 
associated with sea surface heat exchanges. Time-longitude 
maps of departures from the annual mean (denoted with a prime) 
G6feuhe zonal (u') and meridional (v') flow, vertical velocity 
(w'), and temperature (T') were developed to show any zonal 
propagation of wavelike disturbances. Figures 5 and 6 are 


Such maps of u' and v' respectively at level 5 across 12 2S 
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over the last two years of the model run. These figures in- 
dicate westward propagating annual signals which emanate from 
near or at the eastern boundary with the zonal component hav- 
ing a slightly faster phase speed and a higher magnitude. 
ieee / 1S a Similar depiction of the disturbance vertical 
velocity, w'. The major features of the open ocean vertical 
velocity field in Figure 7 develop near the eastern boundary 
and seem to propagate westward as far as the mid ocean region 
at which point further westward propagation with increasing 
time is not observed. White (1977) observed that due to the 
inverse latitudinal qoneaenee in zonal phase speed, baro- 
clinic long waves in the eastern ocean are directed westward 
while in the western ocean they appear to be directed more 
northward due to refraction effects. Regardless, the overall 
pattern observed in Figure 7 does not correlate with those 
of the zonal and meridional flow components shown in Figures 
o and 6 respectively. Vertical velocity is generally related 
to the divergent part of the flow; the propagating features 
See reures 9 and 6 are most likely due to the rotational part 
feet ow and as Such have little impact on the vertical 
memeetry. If 1S noted that disturbances in Figures 5 through 
7 appear to intensify along the western boundary of the model 
ocean. This is to be expected Since westward propagating 
features tend to congregate at or near western boundaries 
in models such as that used in this study. 

Beasineeoms salsO, SHOWS Local disturbances Of great inten—- 
sity along the eastern boundary. However, these disturbances 


do not appear to propagate westward to any significant extent. 
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They also do not appear to be related to the longshore (merid- 
1onal) component of wind stress nor to local Ekman pumping. 
In order to investigate any meridional propagation of these 
boundary features, a time-latitude map of w' along the eastern 
boundary was prepared and is shown in Figure 8. Wyrtki (1975) 
has surmised that the annual baroclinic response at an eastern 
boundary away from the equator iS in part composed of a per- 
turbation which propagates into the region from the equatorial 
region in the form of an internal Kelvin wave. Those coastal 
disturbances at 12.2°N in Figure 8 have a computed (northward) 
phase speed of 11.7 cm/s and the major features have an annual 
period. It is interesting to note that Enfield and Allen 
(1978) have recently observed spatially coherent disturbances 
in monthly mean anomalies of temperature and sea level along 
meme acific coast of North and South America. In the tropical 
latitudes, these observed fluctuations are also uncorrelated 
with local forcing mechanisms and propagate poleward at phase 
speeds comparable to those obtained in the model (Allen, 
personal communication). Because of the coarse grid in the 
present model however, no further comparison between model 
results and coastal observations 1S justified. 

Another important result which can be seen from Figures 
9 through 8, and others below, is that the two consecutive 
model years are virtually identical. This means that the 
model 1S in a statistical equilibrium with the forcing and 


that the oscillations being examined are not transitory but 


are permanent features of the model circulation. 








Figure 9 is a time-longitude map of temperature dis- 
turbances which is very consistent with the vertical velocity 
pattern given in Figure 7. A two month phase lag is apparent 
with generally warm water sinking during the first half of 
the year and cold water rising during the second half - 
clearly and indirect, thermally forced disturbance. Calcula- 
tions based on the annual positive disturbance in the east- 
central ocean in Figure 9 revealed a zonal wavelength of 
about 10,000 kilometers and a zonal phase speed of approxi- 
mately -42 cm/s. These results are in qualitative agreement 
with White (1977) who theorized a wavelength of 10,000 kilo- 
meters and a phase speed of -32 cm/s and observed a wavelength 
of 16,000 kilometers and a phase speed of -46 cm/s at aN 

The Ekman pumping velocity was calculated from the 
following formula as described in detail by DeWitt and 


Leetmaa (1978) 


A map of the disturbance Ekman pumping velocity is presented 
mieeicure 10. This pattern corresponds with the gross fea- 
tures of the vertical velocity field of Figure 7 although no 
propagation is noted in the eastern ocean. This is in con- 
meast with the results of both White (1977) and DeWitt and 
Leetmaa (1978) who noticed a tendency for the Ekman pumping 


disturbances to propagate westward. Magnitude differences 


between corresponding features in Figures 7 and 10 are perhaps 








due to the fact that Figure 10 represents the depth of the 
bottom of the Ekman layer which need not be coincident with 
level 5 of the model. Additionally, different processes 
are involved in vertical velocity and Ekman pumping. The 
theoretical Ekman depth was calculated according to the 


hormu la 


(DeWitt and Leetmaa, 1978) where Kn is the vertical eddy vis- 
weer y coefficient in the model which is 10 om“/s, AMG ds 
the Coriolis parameter. At 12 SON. this depth was calculated 
to be approximately 20 meters. 

In order to describe the vertical structure of the 
memual wave disturbances, vertical cross sections (x,z) of 
the disturbance quantities were prepared across 122 on 
August 3 of the last model year. This date was selected so 
as to coincide with the large positive amplitude disturbance 
noted on the time-longitude section of Ekman pumping in 
Meemwre 10. Figure 11, the vertical cross section of u* on 
August 3, may be compared with the time-longitude map of u" 
in Figure 5. The positive wavelike disturbance near the 
eastern boundary in Figure 5 extends from the surface to the 
bottom with greatest magnitude near the surface, as is the 
case with the positive disturbance in the western ocean. The 
negative disturbance in the mid-ocean extends from the surface 
to about 1500 meters with greatest magnitude at about 50 meters. 
The result that the largest magnitudes are not at the surface 


moreire i] may be explained by the dependence of the zonal 
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flow on temperature and pressure gradient structure and not 
SO much on the surface winds. Generally, then, the Zonal 
disturbances extend quite deep at their formative stage and 
do not Significantly deepen or shoal as they propagate west- 
ward. 

Figure 12 shows the corresponding cross section of 
the meridional velocity perturbation which may be compared 
With the time-longitude section of Figure 6. The positive 
disturbance in the mid-western ocean in Figure 6 on August 3 
extends vertically to within about 950 meters of the surface. 
The negative disturbance in the eastern ocean extends from 
about 1500 meters to the surface where it expands to cover 
the entire model domain across One The negative distur- 
bance is of greatest magnitude near the surface at all longi- 
tudes while the positive disturbance is very small whenever 
it appears. Apparently the meridional flow component is 
more closely related to the surface winds than is the zonal 
eemvonent. 

Promre 15 depicts the vertical=zZonal cross section 
Smevyvertical velocity disturbance across 12.2°N on AUCUSE ok 
The positive disturbance in Figure 7 encompasses most of the 
meade! domain except for the surface layer and the extreme 
eastern boundary. Largest magnitudes are centered at depths 
Of about 1000 meters. The eastern boundary negative distur- 
Dance extends from about 200 meters to the surface. 

Corresponding temperature disturbances are shown in 
Figure 14. Warm conditions are noted in all areas on August 


3 except the eastern region at levels below about 75 meters. 
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MiVschHiwlked (penas noted that lf wind forcing is as— 
sumed to be over an annual period then there is no reason to 
isolate the problem of vertical propagation from that of 
horizontal. He contends that internal gravity waves will 
determine the vertical modal structure and that in the 
presence of a non-dissipative and flat bottom energy may 
travel from the surface to the bottom and then be reflected 
back to the surface. Any vertical propagation of wavelike 
disturbances should be apparent in time-depth maps of the 
perturbation quantities. These depictions were developed 
and are shown in Figures 15 through 18. These maps were made 
at an arbitrarily selected interior ocean longitude of 162°W. 
No vertical propagation iS seen in the zonal and meridional 
velocity components, Figures 15 and 16 respectively. Again 
it can be seen that u' penetrates deeper than v' and is 
therefore most likely due to the north-south meandering of 
the major equatorial current systems which are primarily 
menal. Figure 17 reveals that the vertical velocity distur- 
bances also do not appear to propagate vertically but they 
attain their greatest magnitudes at fairly substantial 
depths. These deep centers apparently are the cause of the 
corresponding deep temperature disturbances (Figure 18) 
feeen lag the vertical velocity by about 2 months. The 
temperature fluctuations above 150 meters are due to the 
seasonal cycle of surface heat flux modified slightly by a 


fluctuating meridional temperature advection in the upper 30 


meters. Again notice how the two model years are nearly 
feentical — thus in equilibrium with the forcing. 
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The model results at 12.2°N are summarized as follows. 
The dominant thermal response at the depth of the main thermo- 
cline occurs in phase across the entire longitudinal extent 
of the ocean and is attributed to large scale vertical advec- 
tion produced by wind-induced Ekman pumping. This is in 
agreement with the theoretical analysis and observational 
results of Meyers (1978). In addition, there iS a pronounced 
fluctuation in the zonal current in the open ocean which 
appears to be associated with a north-south seasonal meander- 
ing of the main Zonal current systems. There also exists 
poleward propagating disturbances confined to the immediate 
Vicinity of the eastern boundary which do not appear to be 
related to the local forcing. Further study is required to 
determine whether these disturbances are in any way related 
to the sea level disturbances observed along the North 
American coast by Enfield and Allen (1978). 

b. Results at 6.1°N 

Variations in the depth of the 14°C isotherm which 
propagated westward at nearly the phase speed of free non- 
dispersive Rossby waves were observed at 6°ON by Meyers (1978). 
He theorized that the waves were generated by the Zonal and 
Seasonal variability in Ekman layer divergence in the 
eastern Pacific and that they propagated nearly freely in 
the western Pacific where he found the forcing was weak. The 


resulting wave lagged the Ekman pumping velocity by slightly 


more than one month at the position of maximum forcing. 








Figures 19 and 20 are time-longitude maps of zonal 
and meridional flow disturbances, u' and v' respectively, 
across 6.1°N at 162.5 m (model level 5) during the same two 
model years analyzed above at 12.2°N. These figures indicate 
rapidly westward propagating annual Signals which appear to 
emanate from the eastern boundary. The zonal component in 
Figure 19 appears to increase gradually in magnitude to a 
maximum in the region slightly west of the mid-ocean and then 
decreases from there to the western boundary. The zonal dis- 
turbances propagate faster than the meridional and faster 
at 6.1°N than at 12.2°N (Figure 5). Figure 21 is a time- 
longitude map of vertical velocity disturbance at the same 
latitude and depth. As was noted for Figure 7 at 12.2°N, 
local disturbances of great magnitude are observed along the 
eastern boundary which do not appear to propagate westward. 
These northerly propagating Kelvin-like disturbances (it can 
Memseen from Figures 19 through 22 that v', T', and w' are 
involved in the disturbance but not u'‘') are depicted in 
Figure 8. Near the eastern boundary in Figure 21, wavelike 
perturbations form and propagate completely across the model 
domain, in contrast with those at 12.2°N (Figure 7) which 
mepear to be nearly in phase across the ocean. The distur- 
bances' amplitudes slowly vary with westward propagation. 

The temperature disturbance at 6.1°N (Figure 22) 
mee essentially the same pattern as the vertical velocity 
disturbance but with a two month phase lag which indicates 
miestemperature is responding to the vertical velocity. Cal- 


culations revealed a phase speed of -32 cm/s, a zonal 
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wavelength of about 9900 kilometers, and an annual period. 
Meyers (1978) observed a corresponding phase speed of -56 
cm/s with a wavelength of 17,000 kilometers at 6 N. 

A time-longitude map of Ekman pumping velocity distur- 
bance across 6.1°N is shown in Figure 23. The theoretical 
Ekman depth as calculated above for 12.2°N was found to be 
approximately 28 meters at 6.1°N. The Patterns in Fieure 23 
correlates well with the actual vertical velocity (Figure 21) 
which in turn correlates with the temperature fields except 
for the region of very slow westward propagation in the mid- 
ocean area. Propagation is very rapid in both the eastern 
and western thirds with greatest amplitudes observed at the 
eastern boundary (125°W). The latter observation iS consis- 
tent with Meyers' (1978) finding that the maximum displace- 
ment occurs at 130°W. The Secon OmesEGEonarecyclonle wind 
stress curl observed at this latitude in Figure 3 was at 
and slightly west of 140°W. 

The zonal-depth structure of the disturbance patterns 
across 6.1°N on August 3 of the last model year is depicted 
in Figures 24 through 27. Figure 24 shows that the positive 
ieersturbance of Figure 19 1S very broad at great depths, 
narrow at mid levels, and widens again in the surface layers 
With slightly higher magnitudes observed at the surface. The 
corresponding negative disturbance 1S widest and has its 
greatest magnitude at about 75 meters. As was discussed pre- 
viously with regard to Figure ll, there 1s a deep zonal 
Suerent fluctuation which seems to be related to the north- 


south displacement of the main current systems which are 
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Menai. the meridional velocity disturbance across on eas 
presented in Figure 25. A relatively weak flow is predomi- 
nant across the model domain except in the surface layers 
where, aS in Figure 12, there is a strong response due to 

the dependence on surface winds. Positive disturbances shown 
in Figure 25 are of much smaller magnitude. 

The zonal-depth map of w' in Figure 26 shows that 
with the exception of a narrow semi-vertical band in the 
western ocean and the Kelvin-like disturbances on the eastern 
boundary, the model domain is predominantly occupied by sink- 
ing (w' < 0) disturbances at this time. The shallow Kelvin 
disturbance extends deeper at 6.1°N than it did at 12.2°N 
(Figure 13) and a phase change is apparent between the two 
latitudes. Figure 27 shows that the model domain is colder 
than normal except for the eastern ocean surface layer and 
a mid-level warm disturbance in the western ocean. 

Time-depth maps of u' and v' at 162°W are shown in 
Figures 28 and 29. Disturbances of greatest magnitudes occur 
at the surface with those of the zonal component extending 
Quite deep. On comparing these figures with their counter- 
parts at 22 Ne Heres io and wee the ditierences in dis= 
turbance magnitudes between the two latitudes 1S obvious; the 
deep response is a little stronger at 6.1°N but the surface 
Ekman flow is nearly an order of magnitude higher at 6.1 N. 
Plots of the disturbances of the zonal and meridional com- 
ponents of wind stress as functions of time at 162°W are 
shown in Figure 30. The large magnitude surface layer distur- 


mamees noted in Figures 28 and 29 are in phase with these 
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Stress disturbances indicating that the wind stress is most 


likely the dominant forcing process of ocean (Ekman) flow 
at the surface and immediately below. The time-depth map 
of vertical velocity disturbance across 6.1°N as presented in 
f@foure 31 iS Similar to that of 12.2°N (Figure 17) in that 
ereatest magnitudes are reached at deep levels. However, in 
Sonerast to 1B. there is a phase shift at intermediate 
levels at 6.1°N. The phase shift probably represents a 
transition region between the upper 500 meters or so, which 
1s responding to the wind, and the remainder of the ocean 
which is reacting out of phase. Figure 32 shows the strong 
response of the temperature field in the upper 500 meters to 
the vertical motion with a two month lag. At the same time, 
Surface disturbances due to seasonal wind changes occur but 
they are confined to the upper 50 meters. At 6.19N v' advects 
cold water northward from the equatorial surface region during 
the first half of the year and warm water southward during the 
second half of the year. 

The model results at 6.1°N are summarized as follows. 
The dominant thermal response in the vicinity of the main 
thermocline (level 5, or 162.5 meters in the model) is a west- 
Parad Propagating disturbance (Figure 22) produced by vertical 
aagvection (w’ in Figure 21) induced by Ekman pumping (Figure 
23). Figure 27 shows that the temperature response seen at 
162.5 meters involves the entire thermocline and is actually 
Strongest at about 100 meters. In addition to the propagating 
thermal response in the main thermocline there is an equally 


pronounced annual cycle of surface temperature (Figure 32) 
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which is produced primarily by meridional advection of mean 
temperature (colder to the south) by surface Ekman currents 


(Figure 29). 
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IV. CONCLUSIONS 


Wavelike features Simulated in the tropical regions of 
an ocean circulation model which are seen as disturbances 
from the annual mean of vertical velocity and temperature ~ 
have been described. These features are considered to be 
permanent oscillations of the model circulation and result 
from variability in wind stress forcing. At 12. 2°N, a 
thermal response to large scale vertical advection produced 
by wind-induced Ekman pumping or surface layer divergence 
occurs practically in phase across the entire longitudinal 
extent of the model ocean. Underlying features of the 
annual Signal do appear to propagate westward. These results 
are in agreement with the theory and observations of Meyers 
(1978). Evidence also suggests a zonal fluctuation in the 
interior ocean current systems. Along the eastern boundary 
poleward propagating perturbations of the vertical velocity 
field exist which may be related to sea level disturbances 
which have been observed by Enfield and Allen (1978). 

At BALIN thermocline displacements propagate westward 
also aS a result of vertical advection induced by Ekman pump- 
ing. Disturbances at this latitude are Significantly higher 
in magnitude than those observed at 12.2°N. AGG LOna ley , 


pronounced sea surface temperature disturbances are apparently 


produced by meridional advection of mean temperature by 








Sum ace EKmMan Currents. Eastern boundary perturbations also 
exist and extend deeper at 6.1°N than at the higher latitude. 

Results and observations in this study were generally 
consistent with accepted theory. It would be beneficial to 
consider latitudes other than 6.1°N and 12.2°N in future 
studies in order to more completely describe the tropical 
wavelike disturbances which were noted. The equatorial 
boundary in the current model may, however, cause unreason- 
able results at lower latitudes than 6.1°N. The model is 
presently not suitable for use for further study of the 
eastern boundary disturbances because of grid size limita— 
tions. However, a careful study of these features uSing a 
higher resolution model having a realistic coastline con- 
meetracion may be valuable in light of the recent findings 
of Enfield and Allen (1978). 

The present model has little or no synoptic "noise". 


ies recognized that this feature may be altered if finer 


resolution or synoptic wind forcing were introduced. 
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